Anderson localization in Bose-Einstein condensates 
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The understanding of disordered quantum systems is still far from being complete, despite many 
decades of research on a variety of physical systems. In this review we discuss how Bose-Einstein 
condensates of ultracold atoms in disordered potentials have opened a new window for studying 
fundamental phenomena related to disorder. In particular, we point our attention to recent exper- 
imental studies on Anderson localization and on the interplay of disorder and weak interactions. 
These realize a very promising starting point for a deeper understanding of the complex behaviour 
of interacting, disordered systems. 

PACS numbers: 



I. INTRODUCTION 

Disorder is ubiquitous in nature and often, even if it 
is only a weak perturbation, it tends to strongly affect 
the properties of many physical systems. The most cel- 
ebrated effect of disorder is probably Anderson localiza- 
tion, i.e the localization of individual particles or waves 
in a disordered energetic landscape. This fundamental ef- 
fect of disorder has been explored in depth already more 
than 50 years ago motivated by the obvious interest 
in the dynamics of electrons in crystals. There imper- 
fections are indeed almost unavoidable, and can lead to 
a dramatic change of the conduction properties of real 
materials Q. 

Anderson localization can be understood as the effect 
of multiple reflections of a plane wave by random scatter- 
ers or random potential barriers. Since the amplitudes 
or the phases are not all identical, as they would be in 
a periodic potential, they tend to produce a suppression 
of the propagation amplitude of the wave at sufficiently 
large distances, even if the reflected amplitude is exceed- 
ingly small. In this case, it can be easily shown that 
at large distances the wavefunction must then decay as 
\<f)(x)\ oc exp where the parameter £ is the localiza- 

tion length [3[ . Such phenomenon appears already when 
the potential barriers are much smaller than the kinetic 
energy associated to the plane wave, a region of parame- 
ters where classically one would expect the disorder to 
produce just a weak perturbation of the propagation. 
Localization can occur not only for plane waves in free 
space, but also for Bloch waves, i.e. for particles or waves 
moving in a lattice potential. It is well known that in this 
case reflections at the potential barriers are already very 
strong even in absence of disorder, but extended Bloch 
states can however exist, provided that the energy of the 
individual lattice sites are equal or shuffled in a periodic 
way. In 1958, P.W. Anderson demonstrated that if the 
on-site energies are instead shuffled in a random way, the 
reflections interfere destructively and the eigenfunctions 
of the system become exponentially localized [1]. 

Much theoretical work has been done about Ander- 
son localization in the past fifty years with important re- 



sults, for example, in assessing how localization behaves 
in finite size systems (the so-called scaling theory of lo- 
calization |4|). Several decades of experimental study 
of condensed-matter systems have on the one hand es- 
tablished the relevance of disorder to many phenomena, 
with recent highlights such as the quantum Hall effect 
5| or the exceptional conduction properties of graphene 
6]. On the other hand, they have not provided con- 
clusive results about Anderson localization itself. This 
is mainly because in real materials there are ingredients 
that go beyond the approximation of a single particle 
in a stationary disordered potential 0, 0], which make 
the observation of Anderson localization difficult. There 
one has indeed to consider also excitations of the poten- 
tial (phonons) and interactions between particles, which 
are typically strong in electronic systems. These effects 
might weaken or enhance localization of particles, de- 
pending the details of the system. For example, it is 
well known that a repulsive elastic interaction between 
electrons can enhance localization [?| , while the same re- 
pulsion between bosons tends to screen disorder and to 
weaken Anderson localization ||. On the other hand, 
it is believed that inelastic processes, such as electron- 
phonon [2| or electron-electron [9j scattering can provide 
a hopping energy that is sufficient to overcome localiza- 
tion. 

To describe disordered condensed-matter systems one 
then needs to fully understand the interplay of disorder 
and interactions. This has however turned out to be an 
extremely difficult task not only in experiments, but also 
in theory pj . There are currently several open questions 
related to this interplay, not just for electrons but even 
in the conceptually simpler case of bosonic systems. For 
example, there is still discussion in theory on how a weak 
interaction can affect Anderson localization [lol - [l7j . or 
on the quantum phases that arise due to disorder and 
interactions in the regime of strong correlations [TJ EH - 

From the experimental point of view, various engi- 
neered systems have recently appeared where transport 
of particles or waves can be strongly affected by control- 
lable disorder as in the case of electrons in solids. Since 
most of the parameters of these systems are highly con- 
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trollable, they can be employed to study in depth the 
basic aspects of the physics of disorder and interactions. 

The first systems in which Anderson localization has 
been unambiguously observed (for a detailed review, 
see for example ref.[23j|) were random photonic media, 
where electromagnetic waves get localized by randomly 
distributed scatterers HH- Electromagnetic waves, 
as well as sound waves [28l430l |. do not suffer from the 
large nonlinearities of electronic systems, and the single- 
particle character of Anderson localization can be clearly 
seen. The main obstacle to the observation of localiza- 
tion in these systems has instead been the presence of ab- 
sorption. More recently, various groups have been able to 
realize also photonic lattices with disorder, where various 
models of disorder can be tested [3ll - l33j . These systems 
allow also to study the effects of interactions, thanks to 
presence of tunable Kerr nonlinearities. However, these 
nonlinearities arc rather weak, and so far only the effect 
of attractive interactions has been clearly identified in 
experiments. 

Localization phenomena have also been studied in the 
context of the quantum kicked rotor, a system that fea- 
tures Anderson localization in momentum space. After 
that seminal theoretical studies have recognized that An- 
derson localization has its clear signature in quantum 
transport processes in periodically driven systems (33 - 
[37j various experiments have been performed, mainly by 
studying the microwave ionization of atoms [H, HI] , or 
ultracold atoms in free expansion interacting with pulsed 
periodic potentials H3, E| ■ The analogous of Anderson 
localization in ID has been studied in depth in those sys- 
tems, and recent studies have also unambiguously iden- 
tified an Anderson transition in 3D [41]. Ongoing exper- 
iments are providing a new i nsig ht on the properties of 
the localization at criticality |42j ]. 

The latest systems that have appeared on the scene 
of disorder studies are ultracold quantum gases of atoms 
in traps, which are interesting in view of the unprece- 
dented possibility of playing with both quantum statis- 
tics and of controlling the most relevant physical param- 
eters, such as temperature, density, dimensionality and, 
very importantly, interactions [HI, [44| . In addition, us- 
ing laser beams one can create potentials of almost ar- 
bitrary shape, ranging from perfectly sinusoidal lattices 
to random speckle potentials [45j]. The atomic systems 
can therefore be exposed in a controlled way to various 
kinds of disordered potentials, thus realizing several of 
the most interesting disordered models studied in the- 
ory. This combination of full control over disorder and 
interaction in a quantum system is offering new possibil- 
ities for finding a solution to some of the open questions 
related to the physics of disorder. 

After few years of research on disorder with quantum 
gases, promising results have already been obtained in 
the context of Anderson localization (for a review, see 
also refs. [46-48]). Two experimental groups have re- 
ported the observation of localization of a non-interacting 
Bose-Einstein condensate in two different kinds of disor- 



dered potentials [HI, [5(| , and other studies on the inter- 
play of disorder and interactions are on the way [SJ [HJ . 
While these studies have not solved all the open questions 
about Anderson localization, they have clearly shown 
that a number of them can be addressed by employing 
atomic quantum gases. The exploration in both exper- 
iments and theory is however not limited to the weak 
interaction regime that is relevant for Anderson localiza- 
tion. For example, from the experimental point of view, 
notable progresses have been reported also in the under- 
standing of the physics of strongly correlated, disordered 
systems [53ll54|. 

In this review we discuss how ultracold atomic quan- 
tum gases are being employed to investigate the physics 
of Anderson localization. The focus will be on the re- 
cently performed experiments on the dynamics in real 
space and in momentum space of Bose-Einstein conden- 
sates in correlated random potentials and in quasiperi- 
odic lattices. We discuss both regimes of Anderson lo- 
calization i.e. negligible atom-atom interactions and that 
in which there is an interplay of disorder and a weak in- 
teraction. We conclude by pointing out the similarities 
and complementarities of atomic and photonic systems 
in the study of Anderson localization, and by outlining 
the interesting research directions for the near future. 



II. DISORDER MODELS AND ANDERSON 
LOCALIZATION 

As already outlined, at the single-particle level the 
term "Anderson localization" is broadly employed to 
name the exponential decay of the eigenfunctions of a 
quantum system that arises in presence of disorder. This 
phenomenon has of course a strong effect on the trans- 
port properties of macroscopic systems, e.g. electrons 
in a crystal or light in a material slab, where the An- 
derson localization is seen as a dramatic change of the 
conductance or of the diffusion constant. The details of 
how such localization occurs in a disordered system is of 
course strongly dependent on the type of disorder and on 
the energy spectrum that the system would have with- 
out disorder. Many models of disorder that are relevant 
to real systems have been studied in the past; their list 
is too long to be even summarized here, but there are 
extensive works on the subject, such as the theoretical 
review presented in Ref. [551] V 

A disordered system that is often considered in theory 
consists of a free particle experiencing a series of uncor- 
related random potential barriers. The barrier heights 
are supposed to have a continuous distribution, with no 
lower bounds. This corresponds to have a <5-like correla- 
tion function for the potential, g(z) = (V(x)V(x + z)) = 
A 2 5(z). Throughout this work the quantity A will repre- 
sent the mean potential energy shift associated to disor- 
der, or disorder strength. For this special case one finds 
that in an infinite ID system all the eigenstates are al- 
ways localized, even for exceedingly small values of the 
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average barrier height. This is a consequence of the fact 
that in ID each barrier causes a complete reversal of the 
motion of a part of the incident wave. The localization 
length £(E) is however approximately linear in the ki- 
netic energy, as one would intuitively expect since the ef- 
fect of the quantum reflections decreases as E is increased 
much above the potential barriers. In ID it is possible 
to establish a series of fundamental relations that help in 
understanding the properties of the disordered systems. 
For example, a relation between the localization length 
and the density of states p{E) that can be useful in many 
instances has been demonstrated by Thouless [56j 
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t(E)- 1 = J \n\E - e\dp{e) . 



(1) 



Also, in the limit of large energy compared to the disor- 
der, E > A, another relation between £ and the correla- 
tion function can be derived [551 
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For the disordered potential above this gives £(E) oc 
E/A 2 . Note that in a random potential the case E < A 
must be treated separately, since there is also the possibil- 
ity of having " classical" bound states in the few deepest 
wellSjJn the so-called Lifshits tail of the energy distribu- 
tion [55l | . 

An analogous disordered potential in 2D gives rise to 
a similar picture, i.e. the states are always localized in 
an infinite system, but the localization length now in- 
creases exponentially with the energy, £(E) oc exp(i? 1 / 2 ). 
This different behaviour can be understood in a random- 
scattering picture as a consequence of the larger number 
of possible paths in 2D, or simply of the fact that a reflec- 
tion process doe not necessarily imply a change of sign 
of the velocity. For the 3D case, theory indicates that 
the states are not necessarily always all localized. For a 
finite disorder strength there exists indeed a phase tran- 
sition between extended and localized states at a specific 
energy E m , the so called mobility-edge. The eigenstates 
of the system are localized only for E < E m , while they 
are all extended for E > E m . For a given system, the 
position of the mobility edge depends on the disorder 
strength. Note that the determination of £(E) is instead 
more complex in higher dimensionality than in ID, and 
often one can only rely on numerical calculations 0, ■ 

As we will outline in the following, the results found 
for this ideal case change when one considers particles 
with different dispersion relations, e.g. particles in a lat- 
tice, or more realistic correlated disordered potentials. In 
the rest of this review we will concentrate our attention 
to the ID case that is more relevant to the experiments 
performed so far with Bose-Einstein condensates. 



A. Random disorder in free space 

A type of disorder that is close to the random potential 
disorder above can be created for Bose-Einstein conden- 
sates by employing laser speckles [H, [llj]. The quan- 
tum gas is initially prepared by evaporative cooling in 
a smooth harmonic potential. Random disorder is then 
added separately by employing the static dipole potential 
of a laser beam that has passed through a diffusive plate 
and has been then focused onto the atomic sample. The 
plate creates a random phase shift along the beam pro- 
file, which is converted into an intensity distribution in 
the far field that has randomly distributed peaks. If the 
laser wavelength is for example smaller than the wave- 
length of the relevant atomic transition, the correspond- 
ing dipole potential is repulsive, and such intensity peaks 
act as potential barriers for the particles, as shown in 
Fig. la. Supposing that the disorder is present only along 
one of the axes of the harmonic potential, one has to deal 
with a ID problem described by the Hamiltonian 



H = -h 2 V 2 /2m 
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V(x). 



(3) 



The relevant information about the disordered potential 
is gained by looking at the energy distribution, which is 
exponential, P(E) cc exp(— E/A) and at its correlation 
function g{z) — (V(x)V(x + z)). Supposing that the 
lens that creates the speckle has a a focal length / and 
a diameter D, a good approximation for the correlation 
function is 



g{z) = A 



2 sin 2 (z/a r ) 
(z/a r ) 2 



(4) 



where A is the mean height of the barriers and a r = 
Xf/irD represents the correlation length of the poten- 
tial. The presence of a non-zero correlation length must 
be considered with some attention, since it changes sub- 
stantially the localization properties of the system in- 
troducing an effective mobility edge, i.e. it allows the 
presence of extended states even in ID (5tl [6(j. As an 
order of magnitude, for a diffraction-limited system one 
has 7roy ~ A, where A «l//m is the wavelength of the 
laser used to create the speckles. How to intuitively un- 
derstand the effect of a finite correlation length? As we 
already mentioned, the localization of a plane wave arises 
from the destructive interference of the reflection from in- 
dividual potential barriers. Now, we can imagine that if a 
plane waves has a wavelength smaller than a r , it will not 
feel the barriers like steep potential edges but more like 
smooth variations of the potential, therefore with a van- 
ishingly small reflection amplitude. If this is the case, the 
plane waves with energy E = h 2 k 2 /2m such that their 
wavevector is k > 1/ay cannot be strongly localized by 
disorder. This point can be seen in a more formal way 
by employing the Wiener-Kintchine theorem, that estab- 
lishes a relation between the correlation function and the 
momentum components of the disordered potential [59j j . 
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Theoretical studies [6l|, l62| have indeed found that for 
a realistic speckle potential there's an effective mobility 
edge at an energy E m = h 2 /2ma 2 , and the appropriate 
form of the localization length within the Born approxi- 
mation is 



8oy EE n 
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(1 - ka r ) 6(1 - ka r ) ' 



(5) 



where Q(x) is the Heaviside function. In comparison to 
the uncorrelated disorder there is then a faster increase of 
the localization length with energy and, at least at first 
order, a strict divergence at ka r = 1. Careful studies 
that go beyond the Born approximation actually show 
that the localization is not completely destroyed above 
the effective mobility edge, but £ increases by at least on 
order of magnitude, depending on the disorder strength. 
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FIG. 1: Cartoon of the spectrum of localized states in the two 
different systems that exhibit Anderson localization realized 
with Bose-Einstein condensates, a) Correlated random dis- 
order in free space, created by laser speckles. The red states 
are exponentially localized, while the blue one is an extended 
state above the mobility edge and the black one is a strongly 
localized state in the Lifshits tail; b) Quasiperiodic lattice for 
A/ J > 2. All the states related to the first energy band of 
the principal lattice (grey area) are exponentially localized. 
Here only the lowest set of states, separated on average by 
d/(/3 — 1), is shown. 



B. Quasiperiodic lattices 

From the discussion above one could imagine that hav- 
ing an indefinite increase of the correlation length would 
lead to a vanishing of Anderson localization. However, 
there are special cases of potentials with non-decaying 
correlation functions that can support localized states, 
the price to pay being the need of a finite disorder 
strength for localization to occur. A paradigmatic exam- 
ple is the quasiperiodic lattice, that consists of a main lat- 
tice perturbed by a weaker lattice with a lattice constant 
that is incommensurate to the first one. In the absence of 
the perturbing lattice, the solutions of the single-particle 
problem are the well known extended Bloch states, that 
are formed by coherent superposition of the Wannier 
states of individual lattice sites. It is easy to imagine 
that if a non-periodic shuffling of the on-site energies is 
added, one might get into a situation where it is no longer 
possible to form extended Bloch states. As we outline in 
the following, theory shows indeed that above a critical 
strength of the secondary lattice the eigenstates of the 
system get exponentially localized as in random disor- 
der. It is interesting to note that an analogous behavior 
appears in quasipcriodically driven systems, where local- 
ization occurs in presence of incommensurate driving fre- 
quencies, see for example ref.[(33| and references therein. 

Crystals with such quasiperiodic lattices, the so-called 
quasicrystals, are well known in solid-state physics [64j j . 
The realization of quasiperiodic lattices for ultracold 
atoms is particularly simple (53|. An individual sinu- 
soidal lattice, the so-called optical lattice, can be cre- 
ated with the dipole potential of a laser beam arranged 
in standing wave configuration. A laser radiation with 
wavevector k gives rise to a lattice with spacing d = n/ki, 
and with a depth U that is proportional to the laser in- 
tensity [45l | . The combination of a deep primary lattice of 
wavevector ki with a secondary one of incommensurate 
wavevector k 2 realizes the quasiperiodic lattice. A sin- 
gle particle in the first energy-band of a potential of this 
kind is actually a realization the well studied Harper (65j 
or Aubry- Andre model 66], which is presented in terms 
of the discrete Hamiltonian 



H = -JT,j(\ w j)( w j+i\ + \ w j+i)( w j\) 
+AJ2 J cos(2nl3j)\w J )(w J \. 



(6) 



The first term describes the hopping of the particle to 
neighbouring sites of the primary lattice, with energy J, 
while the second contains the quasiperiodic shift of the 
on-site energies due to the secondary lattice. The ki- 
netic term of the Hamiltonian has the very well known 
energy spectrum consisting in a whole set of extended 
Bloch states arranged in a band of energy width 4 J. The 
disordering term is parameterized by the wavevector ra- 
tio f3 = k 2 jk\ and the disorder amplitude A, which is 
related to the potential depths of the two lattices by 
A « U 2 /2exp(-P 2 /J Ui ~/E~r), where E r = h 2 k 2 /2m is 
the recoil energy [671 ] . The incommensurability of the 
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two lattices can be maximized by choosing j3 as the ratio 
of two large consecutive elements of the Fibonacci series. 
This leads to the inverse golden ratio, (3 = (y/5 — l)/2. 
Note that in practice, in a finite-size system containing 
N lattice sites it is unnecessary to have f3 irrational. The 
weaker requirement for f3 is instead that it can be ex- 
pressed as a ratio of relatively prime natural numbers 
(3 = P/Q such that P,Q > N. 

It can be analytically demonstrated that the Aubry- 
Andre model admits only extended states below a thresh- 
old disorder strength A = 2 J, and only exponentially lo- 
calized states above [H, [66j]. As already discussed, the 
disordering potential in eqjS] has a non-decaying correla- 
tion function 

g(z) = A 2 /2n cos(2tt/3z) . (7) 

A localization length can however be defined also in this 
case, since the correlation is weak (Hq . A way to identify 
the value of the threshold A t is to note that the Aubry- 
Andrc Hamiltonian is dual, i.e. one can find a trans- 
formation to the momentum space that gives rise to an 
Hamiltonian with an analogous form. The appropriate 
pseudo-momentum basis for the transformation is 

\h) =^exp(i2nPk l j)\wj) , (8) 

3 

and one readily gets the dual Hamiltonian 

h= -£jJY. l {\h)(h +1 \ + \k l+1 )(k l \) + 

+ ¥ A I>os(2^30|fc ; )(fci|- (9) 

This second Hamiltonian has obviously localized states 
where the first one has extended states, and viceversa. 
Since the transition point from extended to localized 
states must be the same for the two, the only possibility 
is that this happens at A/2 J = 1. Numerical solution of 
eq. [6] for f3 = (V5 — l)/2 shows indeed the presence of 
a very sharp transition between extended and localized 
states in an infinite system. Such transitions turns into a 
smoother crossover not only if the system has a finite ex- 
tension, but also if j3 is less irrational, as outlined above 

The localization length of the Aubry- Andre model can 
be derived analytically from e gfT] as £ = d/ln(A/2J), 
and is thus energy independent |66| . Since localization is 
obtained only for A > 2 J, except for a small region of 
disorder strengths A m 2 J, the localization length is com- 
parable of smaller than d, hence quite short. The local- 
ization scenario for the first energy band of a quasiperi- 
odic lattice above threshold, shown in FigfTJ consists in 
N equally-localized states that span an energy 2A + 4 J. 
Their energies keep memory of the quasi-periodicity and 
are therefore arranged in a sort of minibands structure 
[67| . Note that there is not a Lifshits tail in quasiperiodic 
lattices, since all the localized states have a comparable 
overlap with the two neighbouring states. 



It is important to stress again that the Aubry-Andre 
model does not have a mobility edge, but just a criti- 
cal disorder strength above which all states (in the first 
lattice band) are localized. It can be shown that a 
mobility edge appears instead if long-range tunnelling, 
e.g.tunnelling to the next-neighbouring sites, is allowed 
foa . [69[ . This is a regime that can be reached with ultra- 
cold atoms by employing relatively shallow optical lat- 
tices. 



C. Random disorder in a lattice 

It is now instructive to compare the two types of disor- 
dered systems above with the paradigmatic condensed- 
matter disorder model, i.e. the Anderson model [l|. In 
the latter one considers a particle moving in a lattice 
with a random shift of the on-site energies. Simple phys- 
ical systems that realize the Anderson model that have 
been investigated so far are microwave in guides [H, H3| 
and photonic lattices [U [Hj], but this is a configura- 
tion that can be explored also with Bose-Einstein con- 
densates, for example employing a deep optical lattice 
that is perturbed by a weak speckle potential. Experi- 
ments with condensates are actually are underway [54j . 
although not yet in the regime of negligible interactions 
that we are discussing in this review. The ID Anderson 
model is represented by the discrete Hamiltonian 

H=-J ^2(\wj){w j+1 \ + \w j+1 ){wj\) + ^2e j \w j ){w j \, 

3 3 

( 10 ) 

where the first term describes again the hopping and 
the second one contains the on-site energies ej, which 
are supposed to be randomly distributed in the interval 
[— A,+A]. Assuming an uncorrelated random disorder, 
a perturbation theory approach gives in this case a local- 
ization length of the form [7| 

£(E) w 24d(AJ 2 - E 2 )/A 2 , (f f ) 

in the limit of A/J <C 1. The localization length is again 
inversely proportional to A 2 as in free space, but the en- 
ergy dependence is remarkably different, since the max- 
imum of localization is obtained at the two band edges 
E « ±2J, where the energy dependence is quadratic, 
£(E) cx E 2 . Comparison of this result with eq(5] indicates 
that for low energies and the same speckle disorder, the 
localization length in a lattice is smaller than that in free 
space by a factor of the order of E m /AJ, because of the 
lower kinetic energy. It is interesting to note that consid- 
ering the finite correlation length of the speckle potential 
might lead to the appearance of an effective mobility edge 
in ID, in analogy with the free-space case [59j. However, 
in the first band of the lattice the extended states appear 
only if the correlation length is larger than d, and they 
are located close to the band center, in analogy to what 
happens for the extended states above the true mobility 
edge of the 3D Anderson model. 
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D. Anderson localization and weak nonlinearities 

Before discussing the experiments performed with dis- 
ordered Bose-Einstein condensates, we need to briefly 
outline also the effect of a weak nonlinearity on Ander- 
son localization. The fact that nonlinearities such as 
those arising from interactions between the particles can 
change the localization properties of a many-body system 
has been recognized since long, for example in condensed- 
matter systems Q • The basic effect of a nonlinear inter- 
action term in a disordered system is to provide a mean 
of mixing the localized single-particle states to give rise 
to new many-body states. In some cases this results only 
in a small modification of the localization properties of 
the system. There are notable cases where however a 
weak nonlinearity is sufficient to turn the eigenstates of 
the system from localized to fully extended, hence de- 
stroying Anderson localization. The details of how non- 
linearities affect a disordered system depend of course 
from many details such as the quantum statistics of the 
particles (Bose or Fermi) , the single-particle energy spec- 
trum, the temperature of the system and the nature of 
the nonlinearity (attractive or repulsive, short- or long- 
ranged) . 

For atoms, the dominant interaction at ultralow tem- 
peratures is the s-wave contact interaction, which is 
parametrized by the scattering length a. The interac- 
tion energy per particle in a system of N at atoms in a 
state 4>(r) is [zjj 

E int = —aN at f\^v)\ 4 d 3 r. (12) 
m J 

Here a repulsion between atoms corresponds to a posi- 
tive scattering length, i.e. Ei n t > 0, while the opposite 
applies for attraction. In a standard Bose-Einstein con- 
densate such interaction energy can easily be comparable 
or even much larger than the potential energy of the dis- 
order. However, for some atomic species it is possible 
to precisely tune the scattering length from positive to 
negative values across a zero by employing Feshbach res- 
onances (7l| . For example, in a magnetic Feshbach res- 
onance the scattering length can be continuously mod- 
ified simply by applying a homogeneous magnetic field 
that brings different internal states in resonance with the 
atomic state under consideration. 

The case of most relevance for ultracold Bose-Einstein 
condensates is that of a repulsive interaction (i.e a >0), 
since for attractive interactions they are typically not sta- 
ble. Let us then discuss intuitively the effect of an inter- 
action energy that is weak, i.e. is comparable to the the 
potential energy of the disorder, but not much larger than 
the kinetic energy of the system. Let us assume that all 
the single-particle states of the disordered potential are 
localized, as it happens for example in quasiperiodic lat- 
tice above threshold. For zero interaction all the atoms 
will tend to condense in the localized ground state. A 
weak repulsion will instead force the bosons to occupy 



more than one of these states to minimize the system en- 
ergy. These states can further mix between themselves to 
give rise to an extended collection of less-localized states 
if Ei n t is comparable to their energy separation, or even 
to a fully extended system if Ei nt is comparable to the 
full disorder amplitude A. This intuitive picture will be 
further discussed in Section 4, in relation to experiments. 
It is interesting to note that the effect of a weak inter- 
action on a bosonic system is qualitatively different from 
that on a fcrmionic system. In the latter, a repulsion 
tends to keep the particles in their single- particle states, 
which are localized. An attraction on the contrary can 
favour the build-up of extended many-body states that 
span the whole sample, because this allows to establish 
an Ei n t < that lowers the system energy. 



III. EXPERIMENTS ON ANDERSON 
LOCALIZATION WITH BOSE-EINSTEIN 
CONDENSATES 

The idea of using Bose-Einstein condensates to explore 
Anderson localization appeared already several years ago 
[721 [73j , but it took a while before all the required con- 
ditions for its observation in experiments could be met. 
Initially, the focus has been on random disorder in free 
space, created with laser speckles [74Tj77| . These pioneer- 
ing experiments served to identify two different issues. 

The first one was to realize disordered potentials that 
could give the proper range of localization lengths, i.e. 
smaller that the maximum system size, that is typically 
in the range of 1mm, while keeping the disorder strength 
to a minimum to avoid populating just strongly localized 
states in the Lifshits tail of the spectrum. For a speckle 
potential this condition requires a grain size o> that is 
of the order of 1/zm. Since this is the wavelength of the 
laser light used to create speckles, a proper design of 
the optical system is required. While some groups have 
worked in this direction [FJ, [78[ , others have chosen to 
work with quasiperiodic lattices, which naturally provide 
a very small £ [791 ] . 

The second issue was related to the unavoidable pres- 
ence of interaction: all the pioneering studies have ac- 
tually been performed with 87 Rb atoms, which have a 
naturally large repulsion, parametrized by a=100ao, and 
for which the magnetic control of the interaction is ex- 
tremely difficult. The obvious way to solve this problem 
is to reduce theparticles density, either by reducing the 
atom number 17911 or by letting the system expand to a 
large volume |78| . However, the very weak interactions 
achievable in this way can have a non negligible effect 
on the localization, and a real cancellation of the inter- 
action can be achieved only via a Feshbach resonance, 
as we demonstrated for example with 39 K atoms. Actu- 
ally, while some initial experiments were indicating the 
onset of a localization regime, it has not been possible to 
achieve an unambiguous observation of Anderson local- 
ization of atoms in a Bose-Einstein condensate until two 



recent experiments in Paris and Florence [50| . 



A. Localization in a speckle potential 

To observe Anderson localization the Paris team has 
combined a diffraction-limited lens system for the pro- 
duction of speckle potentials with small a r [78j with a 
reduction of the atomic density in a one-dimensional ex- 
pansion of the sample in a waveguide. This has allowed 
a careful investigation of the diffusion of matter- waves in 
a random potential, with a detailed study of the local- 
ization length in ID (49[. 

The Paris experiment employs a Bose-Einstein con- 
densate of 87 Rb atoms, therefore naturally interacting. 
The condensate is initially prepared in a harmonic trap, 
and subjected to speckle disorder with amplitude A/h = 
10 — 100Hz and a r —0.26fim (this corresponds to a mo- 
bility edge at E. m /h = h/2maf « 850Hz). The radii 
of the sample in the trapping potential are about 35/im 
in the axial direction and 3/im in the two transverse di- 
rections. The corresponding E int is as large as /ix220 
Hz, hence definitely non negligible. The idea was then 
to let the condensate suddenly expand out of the trap 
into a one-dimensional optical waveguide aligned along 
the weak axis of the trap and transversely to the direc- 
tion of propagation of the speckles beam. During the 
initial stage of the expansion Ei nt is converted into ki- 
netic energy, the condensate rapidly diffuses out of the 
initial tight spatial distribution, and gradually becomes 
almost non-interacting as its density drops with the ex- 
pansion. At a later stage its wavefunction gets to be a 
superposition of the single-particle localized eigenstates 
of the disordered potential, and the motion therefore is 
confined in space. In particular, the decay of the tails of 
the atomic distribution is determined by the localization 
properties of the largest momentum components, which 
have an energy E Ei nt . To get into the regime where 
trivial localization in deeply bound states cannot be pos- 
sible, care has been taken to have Ei nt > A, while the 
maximum energy has been kept below the effective mo- 
bility edge of eqf5]by having Ei nt < E m . 

Fig. [2] reports the typical experimental observation. 
The axial density distribution of the atomic sample af- 
ter some expansion in the guide is measured by absorp- 
tion imaging. For vanishing disorder strengths A the 
localization length of eqj5] is much larger than the size 
of the observation region L and the expansion appears 
to be ballistic, i.e the rms size of the sample increases 
linearly with time. As A is increased above a critical 
strength such that £,(Ei nt ) < L the size of the sample 
stops to increase after some initial expansion, and the 
atomic density distribution acquires a clear exponential 
decay character in its tails. Note that a fraction of the 
atoms remains trapped in the central region with the typ- 
ical inverted-parabola shape of the original condensate. 
This is presumably a consequence of a trivial localization 
of a fraction of the sample in more deeply bound states in 
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FIG. 2: Geometry of the localization experiment in random 
disorder: a) an interacting Bose-Einstein condensate is ini- 
tially prepared in a harmonic trap; b) the condensate is then 
let free to expand into a ID guide, where it is exposed to the 
speckle potential, c) Experimental observation of localiza- 
tion: for sufficiently large disorder strengths the condensate 
expansion stops, and the tails of its density distribution ac- 
quire an exponential decay behaviour. Figure reprinted with 
permission from 4jj. Copyright 2010 by Macmillan Publisher 
Ltd. 



the Lifshits tail of the spectrum, E < A, something that 
was observed already in previous experiments [7J-|77j. 

The Paris group has then studied the A-dependence 
of the localization length, which is extracted from an ex- 
ponential fit of the tails of the distribution to a function 
of the kind n{x) = uq exp — |x|/£. The experimental ob- 
servations shown in Figj3] feature the expected decrease 
of £ for increasing A and very positively compare to the 
prediction of eqEJ integrated over the density of states 

The regime of energies that are partially above the 
mobility edge E m has also been studied in the experi- 
ment, by increasing the initial interaction energy to have 
ka r > 1. In this case the largest momentum components 
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are not localized and expand indefinitely, producing a 
depletion of the extreme tails of the density distribution. 
The experimental observation is indeed of a faster de- 
cay of the remaining atomic distribution, which as an 
algebraic dependence n(x) — no/\x\ a , with a ps2, as pre- 
dicted by theory |6lj |. 
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FIG. 3: Dependence of the measured localization length on 
the disorder strength. The experimental data (blue dots) 
are in good agreement with the prediction of theory (dashed- 
dotted line). The red region accounts for the uncertainty in 
the determination of E m and Emt in the experiment. Fig- 
ure reprinted with permission from [4{J. Copyright 2010 by 
Macmillan Publisher Ltd. 

A natural question is now whether the reduction of the 
interaction energy that is achieved through the expansion 
is sufficient to neglect completely its effect on the local- 
ization properties. If the measured localization length in 
Fig |3] is used to estimate the particle density, then one 
finds a reduction of the interaction energy by approxi- 
mately a factor of 100 for the largest disorder, or 3000 
for the smallest one. This brings the initial interaction 
energy to about 1Hz or below, hence almost two orders of 
magnitude lower than the mean disorder energy. While 
such a small interaction energy will probably not affect 
the localization properties on the relatively short time 
scale of the present experiment, as we will discuss later 
its effects might be still visible in the long-time dynamics 
of the system. 

B. Localization in a quasiperiodic lattice 

In Florence we have instead employed a Bose-Einstein 
condensate of 39 K atoms, where the interaction can be al- 
most nulled via a broad zero-crossing close to a Feshbach 
resonance (8l| - [83l |. in combination with a ID quasiperi- 
odic lattice. With this approach the observation of the 
localization phenomena does not rely on a dynamical ex- 
pansion of the sample, allowing the investigation of both 



dynamical and static properties of the localized conden- 
sate. The study of the spatial and momentum distribu- 
tions has allowed to verify the expected localization prop- 
erties of a quasiperiodic system described by the Aubry- 
Andre model. 

In the experiment an interacting BEC is initially pre- 
pared in a harmonic potential [50]. A tight ID lattice 
of wavevector k± is then raised using a slow double- 
exponential ramp, optimized to keep the system in its 
ground state; the lattice beams provide also a radial 
confinement, thus realizing a guide for the atoms. A 
second weaker lattice with incommensurate wavevector 
&2 is then raised with a similar ramp lasting about 
100 ms, i.e. a time longer than h/J. The scatter- 
ing length is eventually reduced to a very small value 
a fa O.lao, which corresponds to a residual interaction 
energy that is negligible in comparison with the kinetic 
energy (Ei nt < 0.01J). The ratio of the two wavevectors 
is j8 = 1032nm/862nm sa 1.1972. In a first stage, this 
set-up has allowed to study the spatial diffusion of the 
condensate along the guide, similarly to the Paris exper- 
iment. A main difference is however that in this case the 
diffusion is driven only by the kinetic energy associated 
to the trapping potential. The typical diffusion dynamics 
is shown in FigJU a) when A = the motion is ballistic; 
b)when A < 2J the motion is still ballistic, but with a 
reduced speed, (because of the reduced tunnelling energy 
in the minibands of the superlattice) ; c) when A > 2 J 
no motion is discernible, in accordance with the expecta- 
tion of localized states with £ < d, which is in turn much 
smaller than the initial size of the system, L pa 20d. Inter- 
estingly, the possibility of tuning J independently from 
A, allowed also to confirm that such change of diffusion 
dynamics has the proper scaling with A/ J, as expected 
from eqJS] The analysis of the tails of the spatial distri- 
bution along the quasiperiodic lattice has allowed instead 
to detect the appearance of an exponential localization 
of the far tails of the density distribution for A > 2 J. 

In a second stage we have studied the momentum dis- 
tribution p(k) of the stationary state of the system in the 
harmonic trap. This is obtained by suddenly releasing 
the sample from all potentials, allowing a long ballistic 
expansion driven by the initial momentum, and eventu- 
ally taking an image of the final spatial distribution along 
the lattice direction. An inspection of the evolution of 
the momentum distribution for increasing A elucidates 
the peculiar localization mechanism of quasiperiodic lat- 
tices. In absence of disorder one has the usually narrow 
p{k) corresponding to the q ps quasimomentum states 
of an ordered lattice. For weak disorder, with A < 2 J, 
two momentum components appear in the first Brillouin 
zone, at2fc 2 -2fci = 2(/3-l)fci and4fci-2fc 2 = 2{2-0)k x , 
which are the two Fourier components of the cos(27r/3j) 
term in eqEJ As A increases, other smaller components 
arise, until the whole first Brillouin zone is occupied for 
A > 2 J, implying the presence of states localized over a 
distance £ < d, as expected. 

The momentum distribution gives also information 
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FIG. 4: Diffusion of a non-interacting Bose-Einstein conden- 
sate in a quasiperiodic lattice, a) Cartoon of the experimen- 
tal realization. An initially trapped condensate (left) is let 
free to expand in an optical guide along which is arranged 
the quasiperiodic lattice. The size of the condensate after a 
fixed expansion time (right) depends on the strength A of the 
secondary lattice, and progressive localization is observed for 
increasing A. b) Axial size of the condensate after 700 ms of 
expansion, for three different values of the tunneling energy 
J/h. The inset shows the exponential decay of the tails in 
the localized regime: the blue curve is a gaussian fit, the red 
curve an exponential fit. Figure reprinted with permission 
from [13]. Copyright 2008 by Macmillan Publisher Ltd. 



about the quasiperiodic spatial arrangement of the lo- 
calized eigenstates of the system. When reducing the 
size of the system in the localized regime of A > 2J 
by increasing the harmonic confinement, an interference 
with period 8k = 2(/3 — l)fci clearly appears. This is due 
to the interference of localized states belonging to the 
lowest-lying energy miniband that are indeed sepa- 
rated on average by dj (fi — 1) ~ 5. Id. The phase of such 
interference varies randomly during multiple repetitions 
of the experiment, indicating that the localized states 
have indeed not a fixed phase relation, i.e. they are in- 
dependent. Actually, if the number of states is increased 
above 3-4 the interference contrast rapidly vanishes, as 
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FIG. 5: Stationary momentum distribution of a noninteract- 
ing Bose-Einstein condensate in a quasiperiodic lattice, for 
increasing strengths of the disordering lattice. The broad- 
ening of the distribution provides evidence of a progressive 
localization. 

expected if the relative phases are random. The occupa- 
tion of several localized states testifies that the system 
is not in its ground state, despite the preparation proce- 
dure has been designed to minimize the excitations. This 
is probably due to the diverging tunnelling time between 
states at the onset of the localization, a common issue in 
the preparation of quantum systems that feature insulat- 
ing phases [80, Hj) . 

It is interesting to note that, while the the Aubry- 
Andre model has been confined to theory for several 
decades, soon after the experimental realization with 
Bose-Einstein condensates an analogous experiment has 
been performed at the Weizmann Institute with photons 
propagating in a quasiperiodic -photonic lattice [33[ . Here 
a narrow light wavepacket is injected into one side of the 
quasiperiodic lattice, and the evolution of its transverse 
profile is monitored along the propagation direction. The 
localization transition has been detected by studying a 
quantity that if often employed to quantify the spread 
of a disordered system, the so-called participation ratio 
(PR): PR= (£„ |^ n | 2 )7E„ |0„| 4 . A sharp variation of 
the participation ratio has been indeed observed when 
varying the strength of the disordering modulation of 
the photonic lattice, centered at experimental parame- 
ters corresponding to A/ J = 2 in eq. [5] 



IV. EXPERIMENTS ON THE INTERPLAY OF 
DISORDER AND NONLINEARITIES 

As already mentioned, ultracold quantum gases with 
magnetically tunable interaction offer also the possibil- 
ity of studying the effect of nonlinearities on Anderson 
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FIG. 6: Counting the localized states of a Bose-Einstein con- 
densate in a quasiperiodic lattice. Momentum distribution for 
varying sizes of the harmonic trap (r tra p = y/h/mu), corre- 
sponding to about: a) one state; b) two states; c) three states; 
d) seven states. In the latter case the interference vanishes, 
since the states are phase incoherent. 
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FIG. 7: Transverse localization on light in a quasiperiodic 
photonic lattice, a) Geometry of the experiment, b) Evolu- 
tion of the participation ratio with the strength of the lattice 
modulation. Here the parameter A/C is equivalent to the pa- 
rameter A/ J of the atomic case. Figure reprinted with per- 
mission from [331 ] . Copyright 2009 by the American Physical 
Society. 



localization. In the following we briefly review a recent 
experiment with a BEC in quasiperiodic lattices in this 
direction, and we relate the findings to those of experi- 
ments with photonic systems. 

Let us start by repeating in detail the reasoning of 
Section II. D. We can therefore imagine to have an ideal 
Bose-Einstein condensate that is prepared at T = in 
the ground-state of a quasiperiodic lattice with A above 
the localization threshold. The system is therefore in the 
regime of Anderson localization. The addition of a weak 
repulsion can instead drive the system into at least two 
different regimes. If the repulsion is so weak that the 
average interaction energy per particle Ei nt is smaller 
than the characteristic disorder energy A, its main ef- 
fect is just to prevent the occupation of the ground state 
alone. The system energy is indeed minimized by dis- 
tributing the particles in more than one single-particle 
state, so that the interaction energy arising from cqfT"2l 
is reduced. In this case, one can imagine that the vari- 
ous states populated by the condensate keep their single- 
particle character and are therefore still localized. This 
regime is often indicated in condensed-matter theory as 
that of an Anderson glass [lH , since it features the main 
properties expected for a glass. First of all, this is an in- 
sulating regime, since the system is composed of spatially 
localized states. Second, the energy cost for creating exci- 
tations, i.e. for moving particles from one localized state 
to another, vanishes as the system size tends to infinity, 
since the energy to be compensated is just the one due 
to disorder. 

For increasing repulsion one can expect that the num- 
ber of localized states that are populated increases. 
When Ei n t ~ A there will be a high chance that states 
that are neighbouring in space become also neighbouring 
in energy on the scale of E int . In this case groups of 
neighbouring localized states start to mix and give rise 
to less localized states, or eventually to fully extended 
states. One might also describe this derealization as 
the results of a screening of disorder by the interaction 
One can indeed imagine that most of the parti- 
cles are arranged in the single-particles states to provide 
a spatially-varying interaction energy that almost com- 
pletely compensates the disordered potential, while the 
few remaining particles can freely move on top of such 
combined potential. For identical bosons this picture cor- 
responds to a fully delocalized system. A not too weak 
repulsive interaction therefore is expected to keep bosons 
delocalized, while the opposite, i.e. a stronger localiza- 
tion, would happen for an attractive interaction. 

To verify this conjecture, we have studied how the p(k) 
changes when adding a repulsive interaction, by means 
of the same Feshbach resonance that used to cancel the 
natural interaction. The experimental techniques are the 
same of Section III.B, except for the fact that during the 
loading in the lattice the scattering length is no more 
reduced to zero, but kept to some positive value in the 
range a = 1 — 400ao- Fig. [8j shows the evolution of the 
width p{k) in the A — Ei nt plane. 
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FIG. 8: Delocalizing effect of a repulsive interaction on a 
Bose-Einstein condensate in a quasiperiodic lattice, a) Root- 
mean-squared width of the momentum distribution for vary- 
ing A and Em'- the interaction tends to make the system 
more extended (narrower momentum width) or, alternatively, 
to increase the threshold A/ J for localization. The white line 
is the value of Eint that is necessary for derealization as de- 
rived from a screening argument. The correlation between 
neighbouring localized states (b) increases for increasing in- 
teraction, while their phase fluctuation (c) decreases, and fi- 
nally reaches the one measured below the localization thresh- 
old (blue area). Figure reprinted with permission from [52| . 
Copyright 2010 by Macmillan Publisher Ltd. 



The trend of the data in Fig(5] a is clear: a system 
that is localized for vanishing Ei n t (p(k) is broad) tends 
to become more extended as the interaction energy is 
increased at fixed A. Eventually, for not too large A, 
the momentum distribution goes back to that of a con- 
densate below the localization thershold. Alternatively, 
the diagram can be interpreted as a shift to larger A of 
the localization threshold for increasing repulsion. Other 
information can be extracted from the profiles of p(k). 
For example, all the states that are initially localized be- 
long to the first miniband of the spectrum, which has 
a width A' <C A. As soon as E int becomes equal to 
A' (white line in FigJS] a), a restoring of the interfer- 
ence modulation that signals the reestablishing of the 



phase coherence between neighbouring states is observed. 
To quantify the correlation between the states, one can 
evaluate the spatially-averaged correlation function of the 
condensate g{z) — J dxip* (x)ij)(x + z), which is directly 
related to the Fourier Transform of p(k) by the Wiener- 
Kintchine theorem: g(z) — F~ 1 p(k). The natural dis- 
tance for monitoring the evolution of the correlations in 
a quasiperiodic lattice is of course the mean separation 
between the lowest laying states, z — dj(fi — 1). Alterna- 
tively, one can measure the fluctuations of the phase of 
the interference pattern in the momentum distribution, 
by repeating the experiment many times under the same 
conditions [52j • Both these quantities are shown in Figj8] 
b-c for one particular value of A/ J sa 12 above the lo- 
calization threshold. In the Anderson localization regime 
g{z) is almost zero, while the phase fluctuations are con- 
sistent with having totally random phases in the system. 
As Ei n t increases there is a first regime of Ei nt < A' 
in which the two observables do not change apprecia- 
bly, confirming that the only effect of the repulsion is to 
promote the occupation of a larger number of localized 
states. Only when Ei nt > A' the correlations increase 
and the phase fluctuations decrease, eventually reaching 
the background value measured for a Bose-Einstein con- 
densate below the localization threshold. 

These observations are in nice quantitative agreement 
with a theoretical study of the ground state of the 
quasiperiodic lattice 
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and also in qualitative agree- 
ment with the general phase diagrams derived for random 
disorder [l(| HH and quasiperiodic lattices [86[ . 

One important question at this point is about the dy- 
namical properties of the extended system that is recov- 
ered for large repulsions and A > 2 J: would the system 
expand indefinitely if the confinement is removed, as it 
would do in the superfluid regime A < 2 J? This question 
has been long debated in theory [l]| 03, H3, IsH , where 
a variety of disorder models have been explored mainly 
by numerical methods. Recent studies [l3|, [3 HH agree 
in indicating that an interacting, disordered bosonic sys- 
tem in ID should exhibit an indefinite expansion if the 
interaction energy is initially large enough to mix neigh- 
bouring localized states. The expansion however happens 
in a sub-diffusive manner, with the size of the system in- 
creasing with time approximately as L(t) — A(to + t) a , 
where the exponent a is not only smaller than unity, as 
one would expect for a ballistic expansion, but it is even 
smaller than the value 0.5 expected for normal diffusion.. 
This is a consequence of the fact that the states that are 
instantaneously mixed are limited in number and their 
number decreases with time as the interaction energy de- 
creases with the progressive expansion, see for example 
ref.[lj]- The expectation is of a rather small exponent, 
a R3 0.2 that should be rather independent from the in- 
teraction energy and the disorder strength. The sub- 
diffusive character of the dynamics implies a very long 
time scale for observing a strong variation of the size of 
the system, that for our quasiperiodic lattice amounts to 
several seconds [88| . However, experiments are on the 
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way to test the predicted subdiffusive expansion behav- 
ior. One important aspect predicted by theory is that 
the subdiffusion should continue indefinitely if the initial 
interaction energy is large enough to provide derealiza- 
tion; this clearly rises once more a question about the 
very long time dynamics that could be observed also in 
the Paris experiment. 

The possibility of studying in a controlled way the in- 
terplay of disorder and a widely tunable interaction is for 
sure the most interesting novelty introduced by ultracold 
quantum gases on the scenario of disordered systems. A 
similar capability was previously offered only by photonic 
systems. One example is the pro pag ation of light through 
nonlinear atomic media [171 l89l . |90( • Another very inter- 
esting example is the transverse localization of light in 
photonic lattices, where nonlinearities can be introduced 
thanks to the Kerr effect and controlled independently 
from the other system's parameter by choosing both an 
appropriate value of the nonlinear index of refraction of 
the material and the light intensity [3ll - l33j |. The kind 
of nonlinearity that has been explored so far in experi- 
ments is self-focusing, which corresponds to an attractive 
interaction in the atomic case. Also defocusing nonlin- 
earities can be in principle explored 31], although the 
strength of the nonlinear term in photonic systems is 
much smaller than the one for atoms, and the study of 
focusing nonlinearities is highly favoured. In this sense 
photonic lattices allow investigations that are comple- 
mentary to those with Bose-Einstein condensates. An- 
other complementarity of the two kinds of systems stays 
in the initialization of the system: while in quantum gases 
the initial state is typically a broad wavepacket that has 
overlap with many localized eigenstates, the macroscopic 
size of photonic lattices allows also to take as initial state 
almost exactly a single localized mode, or even a single 
lattice site. 

The effect of nonlinearities on the transverse localiza- 
tion of light has been tested so far in both randomly 
disordered ID and 2D lattices that realize the Anderson 
model jU HH and quasiperiodic lattices [33[ . All the ob- 
servations can be interpreted in the frame of the heuristic 
argument of Section 2.D: in the regime where the linear 
waves would be localized, the self-focusing nonlinearity 
mixes localized modes that are close in energy and space, 
providing in general a derealization of the otherwise lo- 
calized modes [H, [HI- If the linear waves are instead 
not yet localized, either because the A/ J parameter of 
the of the quasiperiodic lattice is below threshold or the 
localization length in a random lattice is larger than the 
system size, then the effect of the nonlinearity is to pro- 
mote localization. This is in qualitative agreement with 
the observed increase of the threshold for a BEC with 
repulsive interaction [52| . A very interesting study in a 
random ID lattice has shown that the effect of the "at- 
tractive" nonlinearity is opposite on the lowest localized 
modes of the band with respect to the highest modes. 
In the former case it leads to an increase of the localiza- 
tion 32j, something that can be understood in terms of 



a negative Ei nt that promotes the occupation of the ab- 
solute lowest localized modes and shifts them out of the 
band. The partial derealization that is instead observed 
in the latter case can be justified not only as a shift of the 
initially localized states towards the center of the band, 
where the states are less localized, but also as the result 
of an effective "repulsive" interaction that results from 
the particular phase structure of the lattice modes in the 
top of the band. 
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FIG. 9: Effect of a self-focusing nonlinearity on the transverse 
localization in a quasiperiodic photonic lattice, a) Evolution 
of the participation ratio (PR) with the strength of nonlinear- 
ity, which is determined by the power level of the propagat- 
ing light itself. Below the localization threshold (A/C=2) the 
nonlinearity promotes localization, while above it promotes 
derealization, b) Delocalization of the wavepacket for in- 
creasing nonlinearity at X/C—3. Reprinted with permission 
from j33|]. Copyright 2009 by the American Physical Society. 



V. OUTLOOK 

We have so far discussed how Bose-Einstein conden- 
sates can be employed to provide a direct measurement 
of the localization of weakly interacting systems in one- 
dimensional random and quasiperiodic potentials. The 
two kinds of systems provide complementary informa- 
tion on the mechanism of Anderson localization and on 
its interplay with interactions. The possibility of tuning 
the atom-atom interaction in a controlled way through 
Feshbach resonances is an invaluable tool to evidence 
Anderson localization and to explore in a quantitative 
way the effect of nonlinearities. Obviously, also ultracold 
Fermi gases can be employed for this kind of studies; also 
there the interaction can be magnetically controlled, or 
can even be switched off by using a sample of identi- 
cal fcrmions. Let us now outline some of the disorder- 
related phenomena that have been debated over the last 
few decades, and whose details could be finally disclosed 
in upcoming experiments with quantum gases, in both 
regimes of weak and strong interactions. 

On the weak interaction side, both Bose and Fermi 
gases might allow to extend the study of Anderson local- 
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ization in random potentials to higher dimensionality, by 
means of a rather straightforward extension of the tech- 
niques employed so far. For example, there is particular 
interest in studying the characteristics of the mobility 
edges at the crossover between 2D and 3D in presence 
of a controlled interaction [9l], |92|. This kind of stud- 
ies would be interesting also in relation to the observed 
metal-insulator transition in disordered 2D electron gases 
. First experiments in 2D disordered systems have ac- 
tually been started [93] , and others will soon follow also 
in 3D. A question of general interest is also how much 
superfluidity or superconductivity persists in presence 
of disorder. Theoretical studies in this direction with 
Bose- Einstein condensat es I94B and with Fermi gases in 
the BEC-BCS crossover [95l497| are on the way, and fu- 
ture experiments might provide insight in the physics of 
" dirty" Bose and Fermi liquids that is relevant to various 
condensed-matter systems. 

A particularly important question that is directly re- 
lated to the ongoing experiments is about the nature of 
the many-body excitations of an interacting disordered 
system. A theoretical finding is that a system with an 
extended ground state might have excitations which are 
localized by disorder, although with a £ different from 
the the single-particle ones [98[ . Recent studies show also 
that such "many-body localization" may be destroyed is 
the density of excitations are increased above a threshold 
by increasing the temperature of the system. This should 
lead to a temperature-driven metal-insulator transition 
in both Bose [99( and Fermi |9J systems. 



The other interesting limit is the one of interactions 
that are strong enough to drive a lattice system into the 
regime of strong correlations, a regime that we did not 
discuss at all in this review. Here already in absence of 
disorder the interaction alone is capable of bringing the 
system int o the Mott insulating pha s e, a s observed in 
both Bose [lOOf and Fermi gases [lOll Il02| . Despite ex- 
tensive theoretical studies, there are still open questions 
about the interplay of disorder and interaction in this 
regime, even for the bosonic case. For example, at the 
Mott transition the disorder can in principle have two 
opposite effects, either a destabilization of the Mott in- 
sulator towards the superfluid, or a destabilization of the 
supernuid towards a new insulating phase, the so-called 
Bose glass. In theory there's still debate ab out t he ex- 
act shape of the phase diagram [3413, Il03| . The 
experimental investigation of the phas e dia gram of these 
systems has just started [H, H3, Il04| (see [46| for a re- 
view). Analogous phenomena are of course expected in 
the fermionic case. With the refinement of experimen- 
tal and theoretical techniques, the study of even more 
fragile disorder-related phenomena appears to be within 
reach with ultracold quantum gases, such as quantum 
spin glasses and disordered magnetic phases (4§j . 
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